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The Dissipation of Energy in the Tides in Connection with the 

Acceleration of the Moon's Mean Motion. 

By E. 0. Street, M.A., The University, Liverpool. 
(Communicated by Sir Joseph Larmor, F.R.S. Eeceived February 12, 1917.) 

1. In the following paper an expression for the rate of dissipation of 
kinetic energy in a rotating sea is obtained in the form of an integral of a 
function of the surface-current velocities. This expression has been used to 
find an approximate value for the rate of dissipation in the Irish Sea, for 
which much tidal information is available, and the order of magnitude of 
the result obtained suggests that this direct effect of viscosity may account 
for at least an appreciable part of the earth's secular retardation. 

The water is supposed of uniform density, and at the bed of the ocean to 
be relatively at rest. The motion, which, as usual, is regarded as small, is 
supposed non-turbulent, and simple-harmonic with respect to the time with 
a period of 12 hours. In actual tidal motion eddies are present, so that the 
result here obtained for the rate of dissipation may be regarded as a lower 
limit. 

The question of stability is not considered. This would require a 
knowledge of the slight variations of density. It is frequently found* that, 
at a • given station, the maximum density is not at the bottom or the 
minimum not at the top. 

The ordinary equations of hydrodynamics for a flat rotating sheet of 
water are applied to an elementf of the ocean. Provided that the greatest 
linear dimension of the element is small compared with the radius of the 
earth, the error introduced is unimportant. In any element the bottom is 
supposed flat, and the linear dimensions are further restricted by the 
condition that, in the state of relative equilibrium, the depth, to a first 
approximation, is constant. As the results are obtained without introducing 
lateral boundary conditions, the method may be applied by summation to a 
region of any area. 

2. Consider a plane element rotating with constant angular velocity &> 
about a vertical axis, which is taken as the z co-ordinate axis, the bottom of 
the water being the plane 2 = 0. Let the x- and y-axes rotate in their 

* Many examples of this are given in the Reports on the Hydrographical Investiga- 
tions in the North Sea (North Sea Fisheries Investigation Committee). 

t The expression "element" is used to denote a portion of the ocean of such 
dimensions that the conditions here stated are satisfied within it. 
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plane with the given angular velocity &), and denote by IT, V the velocities 
at time t, along and relative to these axes, of the particle of water which at 
that instant occupies the position (x, y, z). Then, if the conditions stated 
above are satisfied, and if the vertical acceleration of the water is neglected 
in comparison with g, the equations of relative motion may be written* 



3V + 2a>U = -g ~ (£««*) + «-AV 



~N 



h (1) 



dt y dy 

where v is the kinematic coefficient of viscosity, A is the Laplaeian operator 

d 2 /dx 2 +d 2 /dy 2 + d 2 /dz 2 , and 

Sf*-Tl/g (2) 

is the elevation of the disturbed surface above its position of relative 
equilibrium, II being the disturbing potential. 

If II is proportional to e la \ and f is supposed a function of x, y only, let 

(U, V) = (u, v) &°\ (3) 

Then u, v are functions of x, y, z only, and the equations (1) become 



3f 

J dx 



•\ 



iav + 2cau = — g ^f- + vAv 

dy 



Y- (4) 

at r y 



J 



For free motion in the absence of viscosity the elimination of £ between 
equations (4) and the equation of continuity leads to 

Al+ ^r~y u ' v) = ' (5) 

where h is the depth of the water in the element, and Ai is the operator 
d 2 jdx 2 + d 2 /dy 2 . In free motion with v not zero the equations (5) will, of course, 
not be satisfied, but the right-hand sides will be functions of v, u } v which 
vanish with v. It may, therefore, be supposedf that in a free oscillation of 
the sea 

Ai (u, v) and — (it, v) 

are of the same order of magnitude. 

If X2 is the angular velocity of the earth's axial rotation, for an element of 
the sea in latitude X the appropriate value of © is 12 sin X or, say, 0*8 X2 for 
the neighbourhood of England. Also for a semi-diurnal oscillation a = 2X2. 

* Lamb, ' Hydrodynamics,' §§206, 316. 

t The validity of this supposition is verified later (§3). 
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Hence 



v_ a 2 -Am 2 _ . 72 v& _ 9 .g x ]0 _ 12 
a* gh gh 



for a depth of 1 metre, taking v — 0*018 C.G.S. units. 

It appears then that of the terms in u, v on the right-hand sides of 
equations (4) the parts vAi % z>Ai v are absolutely insignificant compared with 
the terms in u, v occurring on the left-hand sides of those equations. Thus 
the equations may be replaced by 

ox CZ J 



Hence 



dy 



y 



(6) 



dz 2 



j 



AL 

,dz 2 



0-4-2(0 w x u , v u \ c. 



V 



v \dx . By/ 1 



a 



S_ 2 er-2©> 

kBz 



T — i 



^;< 



ff/d 



io— cv) =£ _—i — 

v \ox dy 



•lor 



^ 



(7) 



^ 



The conditions to be satisfied by u, v are that there is no velocity at the 
bottom and no shearing stress at the free surface (supposed smooth). Hence 

at z = 0, u = v = 0, 

at z = h, cufdz = Bv/dz = 0, 

The solution of the first of equations (7) is therefore 

«. LV = n (i - cosh i (A =f) v 7 K g ± 2fl, )/ y 3 > ) ( d . + 1 1) r (8) 

o-+2a>\ cosh {h^/[i(cr + 2cD)/p]} J \dx dy) 

(o- + 2cd)/p = 2a 2 , (a-2co)/v = 20®.. (9) 

Then with the values of <r, a), v used previously, 1/a = 12 cm., 1//3 = 35 cm., 
so that if A 'is more than a few metres tanhaA and tanhjSA may be replaced 
by unity without appreciable error. 

With these simplifications equation (8) becomes 



u 



Let 



ii + w 



# 



1— g-a(l + *)^ 



<j-f-2ci> 

and in the same way from the second of equations (J) 

<0 



A 4. ~)f 

8^ s ~9y, ? 



(10) 



%—iV 



a — 2ft> 

Putting a = A, these equations give 






(11) 



^o±^o = 



3 4. 3 W 
a- + 2o> \9# dy/ ~ 



V 



(12) 
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to the degree of approximation already used, where u , v are the surface 
values of u, v, and all the upper or all the lower signs of the ambiguities are 
to be taken. 

Combining equations (9), (10), (11), (12), the equation 

U±iV = (u ±LV )(l— exp{— z x /[i(o-±2G))/p]})e l * t (13) 

is obtained for the complex relative velocities U, V. 

3. It is now possible to verify the assumption that in a sea with viscous 
quality 

&i(u,v) and (u, v) 

are of the same order of magnitude. 

The equation of continuity for free motion is now 

^ = -jo(l + |H (14) 

Integrating both sides of equations (10), (11) with respect to z and neglect- 
ing e~ ah , e~& h as before, 



£<*+">* = ^< 1+A) (I + '|! 



(u—tv)dz = — ™-~(l + B) (■= £5- If, 

<t— 2co K ; \ox • ay J 

where Ax = B/3 = (*-l)/2&. (15) 

Hence from equation (14) 



Therefore 



cr-f-Zco cr — AG), 



l + ^ A +f^B)A 1 r+^g^r=0. (16) 



By equations (15) the coefficient of Aif differs from unity by a small 
quantity proportional to v%. 

From equations (10), (11) it follows that u, v are also solutions of 
equation (16), so that 

Aiu-i - u and A x v-. - — v 

gk gh 

are each nearly equal to — 1 ; and the assumption that they are of the same 
order of magnitude as unity is justified. 

2 E 2 
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It is important to notice that since, by equations (10), (11), 



3 2 / i \ °" 
?r— (u + tv) = — exp 

cb 2 z/ 



a , a 



vd 2 u/dz 2 and vd 2 v/dz 2 do not contain a factor proportional to a power of v. 

4. If (p, q, r) denote the real relative velocities of the particle of water at 
(x, y, z) the rate of dissipation of energy in the element is 

f Jf <E> dx cly dz, 
integrated through the volume considered, where 

*="{ 2s (!) !+s (! + i)'}' 

jub being the coefficient of viscosity. For oceanic tidal waves the important 
terms in <j> are 

-{(DM!)' 

Integrating by parts with regard to z, the rate of dissipation per unit surface 
area is found to be 



C h fd 2 p c 2 q\ 7 



(17) 



the integrated part vanishing at both limits. 

Writing u = % + i%, v = v 2 -\-lv 2 

where %, u 2 , Vi, v 2 are real, there follows from equation (3) 

U = (it! cos at—u 2 sin at) + i (u\ sin at + u 2 cos at). 

(J X) • o 2 

Hence p ■— = (ui cos at — %i 2 sin cr£) -^- (% cos at — u 2 sin cr£), 

of which the mean value over a whole period is 

The mean rate of dissipation of energy per unit surface area is therefore, 
from equation (17), given by 

%h f c 2 Ui d 2 u 2 d 2 Vi d 2 v 2 





F = -^ \A Ui -te +U2 M +Vi w +V2 ws dz 



= —\p 



* r av cV. 7 



where K denotes " the real part of," and u f , v' are the conjugate complexes of 
u, v. Hence 

F = — Ifju | E -l (u + iv) -^ (u' — lv')-\-(u — iv) -^ (u' + w') Vdz. (19) 
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Let uq + ivq = P, uq — ivq = Q, and let P', Q' be the conjugate complexes of 
P, Q. Then from equations (10), (12) 

so that u' — cv' = P'{1 — e-** 1 -'**}. 

Therefore 



j o E {^ + w)^K--w')} 



■A 



= -2PPV f E{[l-e— ( 1+t )^^- ft < 1 -')^}^ 

= — 2PP'aB {t/(l — t)} approximately 
= -PP'«. 

Hence from equation (19) 



VIttI+M 



v 2*/ 

or, since /* = z/p, where p is the density of the water, 

E = y x /(2v){VV\/(a + 2co) + QQ' v /(a-2G>)}. (20) 

Let uo = TJ&*, v = V&+ (21) 

where U, V, 0, </> are real functions of x, y. Then the components of the 
relative surface velocities, being the real parts of 

JJ e i (9 + <ri) ~y e i(<!> + crt) 

are TJ gos (0 + at), V cos (</> + at). (22) 

The hodograph of the motion is thus an ellipse,* whose dimensions and 
orientation are functions of the position of the point on the surface. 
Changing the origin of time, the component velocities may be written 

U cos (at' + yjr), V COS (at' — y}r) 

where 2\}r = — <£. 

If therefore W is the resultant relative surface velocity 

2W 2 = U 2 + V 2 + U 2 cos2(ar + ^) + V 2 cos2(^ , -~^) 

= U a + V 2 + (U 2 + V 2 ) cos 2f cos 2^'-(U 2 --V 2 ) sin 2^ sin 2<rt' 
= U 2 + V 2 + (U 4 + V 4 + 2U 2 V 2 cos 4^)4 cos (2atf + x )> say. 

* For drawings of the ellipses at several stations, see " Hydrographical Observations 
in the North Sea," in *the 'Bulletin Hydrographique,' 1910-11, 1912-13 (Conseil 
International pour FExploration de la Mer). 
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Hence, if W, w are the maximum and minimum values of W 9 

W 2 + <u? = TP + V 2 1 

y. (23) 

Ww = UVsin2^r = VVsm(0-cl>)] 

It is convenient to suppose that W is always positive, w may then be 
negative. 

Now PP' = \u + IVq 1 2 

= (U cos 6- V sin <£) 2 + (U sin + V cos <£) 2 

= U 2 + V 2 - 2 UV sin (<9 --</>). 

Hence PP' = (W~-w) 2 "l 

? . (24) 

and similarly QQ' = (W + w) 2 J 

Thus from equation (20) the mean rate of dissipation of energy per unit 
surface area is given by 

F = hp x /(2v){(W^iay^(a + 2co) + (W+wf^/( ( T--2co)}. (25) 

Actual measurements of current velocities show that w is usually very 
small compared with W, so that it might, without serious error, be put 
equal to zero. It is, however, easy to show that, whatever value w may 
have, the order of magnitude of F is unaltered by neglecting w. 

With the values of a and w used before 

x /(<t + 2g>) = v" (3-612), v /( --2fl>) = v/(0-4Il). 
Hence F is proportional to 

or to W 2 — Ww-\~w 2 . 

Denote this expression hjf(w). Then the maximum and minimum values 
of f{w) are respectively 3W 2 and |W 2 , since w must lie between — W 
and W. Also /(0) = W 2 . Thus f(w) must lie between |/(0) and 3/(0), 
and is therefore always of the same order of magnitude as/(0). 

Writing w = 0, <7 = 212, « = 12 sin \ = 12 cos 0, where X is the latitude 
and the co-latitude of the place considered, equation (25) gives 

F = ipW 2 v / (i/X2){ v / (l + cos^) + v / (l-cos^)}, 
or F = !/> WV(vX2) cos |X. (26) 

The Bate of Dissipation of Energy in the Main Fairway of the Irish Sea at the 

Time of Spring Tides. 

5. The expression obtained in equation (26) for the mean rate of 
dissipation of energy by viscosity per unit surface area has been applied 
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to that portion of the Irish Sea which may be described as the Main 
Fairway. The boundary of the area used has been determined approxi- 
mately, so that all regions in which the range of tide is more than about 
one-tenth of the depth have been excluded. The boundary so obtained 
follows the coastline fairly closely, except in Liverpool Bay and between 
Walney Island and the Isle of Man. The unit of area taken for the 
purpose of summation is bounded by meridians at an interval of 15' and 
parallels at an interval of 10'. On a Mercator chart of this region such a 
unit is approximately square. In each of these units (or occasionally in a 
portion of a unit) a value is assigned to the maximum surface-current 
velocity, and the corresponding rate of dissipation is evaluated. 

The values of these maximum velocities at a large number of stations are 
given in the Admiralty manual, ' The Tides and Tidal Streams of the 
British Islands/ These, supplemented in a few cases by values taken from 
their chart of the Irish Sea, enable a value of W 2 to be interpolated in each 
of the units used. As the observed velocities are given in knots and 
quarter-knots, it is convenient to express W 2 in (| knots) 2 . 

In the few cases in which the boundary of the region precludes the use of 
the whole of a unit, the value of "W 2 has been reduced proportionately. The 

Mean Yalues of W 2 in the Units of Area. 
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letters a, h, c, d following a number in the Table* denote that the N.E., KT.W. 
S.W., S.E. quarters of the unit have not been included. 

The area of the unit depends on its position, and to a sufficient approxi- 
mation is equal to 150 cos X square nautical miles, X being the mean latitude 
of the unit. Without appreciable error the portion of the sea considered 
may be divided into four regions : 52° 1ST.— 53° K, 53° K — 54° K, 54° K — 
55° 1ST., 55° N". — 55° 20' N., within each of which one value of X is used. 

From the Table of values of W 2 the following summary of results has 
been constructed : — 



Mean latitude (X). 


2W 2 . 


2W 2 cos A cos |X. 


Number of units (1ST). 


N cos X. 


52° 30' 

53 30 

54 30 

55 10 


4412 

2766 

1366 

686 


2406 
1469 

705 

347 


35f 

40 
20i 
' 7-*- 


21-76 

23-80 

11-76 

4 *14 



Thus 2 W 2 cos X cos |X = 4927 

SNcosX^ 61-46 

where the summation now extends over the whole of the main fairway of the 
Irish Sea. 

Taking p = 64, i; = r4xl0"" 5 , XI = 7r/(12 x 60 2 ) in foot-pound-second 
units, equation (26) gives as the mean rate of dissipation of energy over the 
whole area considered 

5*1 x 10 9 foot poundals per second. (27) 

The area within which this dissipation takes place is 150 2N" cosX. 

or 9220 square nautical miles. (28) 

6. Neglecting the relatively small vertical velocity, the kinetic energy of 
the disturbed motion is 



\? 



= i P 



{p 2j r(f)dz per unit surface area 



{ (ui cos at — u 2 sin at) 2 + (^i cos at—v 2 sin at) 2 }dz. 



Hence, if T denotes the mean value of this kinetic energy taken over a 
whole period, 



T = \p \ (tii 2 + u 2 2 + vi 2 + v 2 2 )dz 
Jo 

^ ip {(u+iv)(u' — i,v') + (u—tv)(u' + iv')}dz. 



(29) 



* The publication of this Table has been sanctioned by the Hydrographer. 
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Now 





(u + iv) (u f — iv f ) dz 

h 



PP' [ {l-.e-*<i+0*} {l-e-^-^^dz 



= PP' 



<h 







1 — 2e~ az cos a# + e~ 2az } dz 



= PP' (A— — ) to the approximation already used 

= PP'A very nearly. 
Hence from equation (29) 

T = £pA(PP' + QQ'), 

so that from equations (24) 

T = i P h(W 2 + w 2 ). (30) 

A rough approximation to the mean kinetic energy of the portion of the 
Irish Sea already considered may be obtained by supposing k and W 2 + i# 2 
constant, and equal to 250 feet and 5 (knots) 2 respectively over the whole 
9220 square nautical miles. With these assumptions equation (30) gives for 
the mean kinetic energy 

9 x 10 15 foot-poundals. (31) 

Owing to the action of external forces, this mean kinetic energy is nearly 
constant over long periods of time. Suppose that in the absence of these 
forces it could be expressed by an equation of the form 

so that t is the period in which the mean kinetic energy would be reduced 
in the ratio e to 1. Then 

\ dt 

= 1 : 8 x 10 6 seconds, 

or about two hours, using the figures in (27), (31). 

For a non-rotating sea of uniform depth, the approximate formula of 
Hough* is, in the present notation, 

v av 

= 5 x 10 5 seconds, 
using the same values of a, v, h as before. 

* < London Math. Soc. Proe.,' vol. 28, p. 275 (1897). 



9. 
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The agreement between these two values of r is as close as could have 
been expected. 

The Effect of Tidal Friction on the Earth's Axial Rotation. 

7. About three-quarters of the earth's surface is covered by water; so 
that if energy were dissipated throughout all the oceans at the same rate 
per unit surface area as in the Irish Sea, the mean rate of dissipation due 
to tidal friction would be 

6 x 10 13 foot-poundals per second. (32) 

This value is probably several times too high, as in general much smaller 
current velocities are observed in the large expanses of water than in the 
confined regions of the Irish Sea. 

8. It is well known that observations of the moon's position indicate an 
acceleration* in its orbital motion referred to terrestrial time units. This 
acceleration might be explained by assuming that the earth's axial rotation 
has a retardation of about 4' of arc per century per century. 

If T is the kinetic energy of the earth's rotation, and I is its moment of 

inertia, then 

T = JIO 2 

, , , 1 8T 2 3X2 

sothat Ti = sr 

Let —d£lfdt correspond to the retardation of 4' per century per century. 
Then 

41 = ^s^ 1002 • 3652 • 24 • 602rl 

where the second is the unit of time. Taking T = 5 x 10 30 foot-poundals, 
the rate of decay of energy is given by 

— _ = 1-6 x 10 13 foot poundals per second. (33) 

ot 

It will be seen that this rate of dissipation is about one-quarter that given 
in equation (32), No agreement of the two results except in orders of 
magnitude could have been expected. 

9. It is interesting to find the maximum surface-current velocity (assumed 
uniform over all the oceans) which will give rise to the retardation of 4' 
per century per century in the earth's rotation. 

* The actual amount of this irregularity is a matter of controversy, but its order 
of magnitude is of the quantity used here. For references, cf. Sir J. Larmor, 
' M. N. K. Astron. Soc., J vol. 75, p. 211 (1915). It is there shown that if the interaction 
between the earth and the moon is considered, any loss of energy is divided between the 
two bodies in the ratio 14 : 1. 
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If the maximum velocity at every point is W, and if at any latitude the 
area of the land is a third that of the sea, the mean rate of dissipation of 
energy would be, from equation (26), 



r k* 



\it 9 W 2 KVO a ) cos i* cos ^dX = irp W 2 E 2 v /(2i/H) 
where E is the radius of the earth. Hence from equation (33) 

•TTTO X \J X 1U 

or W = 2 feet per second = 1*2 knots. 

My thanks are due to Sir J. Larmor for the interest he has taken in the 
work, and to Dr. Proudman for help and advice on many points which have 
arisen. 



The Complete Photo-electric Emission from the Alloy of Sodium 

and Potassium* 

By William Wilson, Ph.D. 

(Communicated by Prof. 0. W. Richardson, F.E.S. Received February 17, 1917.) 

The experimental work described in the present paper suggested itself to 
the writer in connection with an earlier investigation*)- on the law governing 
the temperature variation of the complete photo-electric emission from a hot 
body, i.e. the photo-electric emission from a body in equilibrium with the 
full (black body) radiation corresponding to its temperature. By making 
use of hypotheses contained in the quantum theory, the writer obtained the 
following expression for the current per unit area 

C = AT(l + 2kT/cj> + 2k 2 T 2 /cj> 2 )e~<f>Iw, 

where <f> is the work done in removing an electron from the hot body, and is 

equal to hv 9 v being the lowest frequency of the radiation capable of 

producing a photo-electric emission, and h being Planck's constant. The 

quantity h is the " gas constant " reckoned for one molecule, and A is a 

quantity independent of T, and characteristic of the substance. As the 

expression inside the brackets in the above formula does not differ 

appreciably from unity, the latter is substantially of the same type as 

Eichardson's equation 

C = ATV-*/**. (1) 

* The expenses of this research were partly defrayed by the aid of a grant from 
the Government Grant Committee of the Eoyal Society, to whom my thanks are due. 
+ W. Wilson, < Annalen der Physik, 5 vol. 42, p. 1154 (1913). 



